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Abstract 

Recent claims of a line in the Fermi-LAT photon spectrum at 130 GeV are suggestive of dark 
matter annihilation in the galactic center and other dark matter-dominated regions. If the Fermi 
feature is indeed due to dark matter annihilation, the best-fit line cross-section, together with the 
lack of any corresponding excess in continuum photons, poses an interesting puzzle for models 
of thermal dark matter: the line cross-section is too large to be generated radiatively from open 
Standard Model annihilation modes, and too small to provide efficient dark matter annihilation in 
the early universe. We discuss two mechanisms to solve this puzzle and illustrate each with a simple 
reference model in which the dominant dark matter annihilation channel is photonic final states. 
The first mechanism we employ is resonant annihilation, which enhances the annihilation cross- 
section during freezeout and allows for a sufficiently large present-day annihilation cross section. 
Second, we consider cascade annihilation, with a hierarchy between p-wave and s-wave processes. 
Both mechanisms require mass near-degeneracies and predict states with masses closely related to 
the dark matter mass; resonant freezeout in addition requires new charged particles at the TeV 
scale. 



1 Introduction 



The existence of dark matter (DM) is one of the strongest pieces of evidence for physics beyond the 
Standard Model. Searches in cosmic rays for evidence of DM annihilation or decays are a cornerstone 
of the experimental effort to detect DM. Monochromatic photon lines, though in most models a 
subdominant signal, provide a particularly clean astrophysical signal [H El [31 [2J [5] . 

Several recent analyses have claimed evidence for a distinct spectral feature in the Fermi-Large 
Area Telescope (LAT) [6j photon spectrum at around 130 GeV [3 02] j in regions near the galactic 
center. Evidence for this feature has also been reported in galactic clusters [11] and in non-associated 
sources [15] . although the latter claim remains contentious [13j [T31 [T5] . While the statistics are limited, 
the morphology of the signal may favor an explanation as annihilating DM [161 \T7\ [TT] . The presence 
of an additional photon line at approximately 111 GeV [101112] is highly suggestive, if also statistically 
limited, and would lend more credence to a particle physics explanation [18] . The Fermi collaboration's 
own search for photon lines uses slightly different search regions and methodology and sets an upper 
limit marginally in conflict with the claimed signal [19] , 

It remains to be established whether the excess is instrumental, astrophysical, or representative 
of an overly optimistic characterization of the systematic uncertainties in the galactic background 
[201 121] . However, standard thermal WIMPs are not capable of explaining the Fermi signal, and it is 
of interest to work out the necessary structure in DM models which could give rise to the 130 GeV 
line. 

Any dark matter model for the Fermi 130 GeV line must account for two interesting facts. First, 
there is no evidence for an excess in the continuum photon spectrum, which strongly constrains DM 
annihilation into the usual SM annihilation channels (//, W + W~, ZZ) [TTJ [55J [23] [52] or indeed 
into any charged final states. Since in most thermal models DM-photon couplings are generated 
radiatively from the DM couplings to charged final states [251 126] . the typical line cross section is 
generically related to the cross-section for annihilation into charged modes X, X' by 

/a\ 2 

(ov)yy ~ {-) (<rv) X xi ■ (!) 

As the fragmentation and decay of the final states X, X< give rise to a continuum photon spectrum 
d^ 1 {E) / dA oc n 2 DM ((Tv) X xi dN^ j 'dE ', if annihilation to charged states is open, the expected line flux 
is smaller than the continuum flux by several orders of magnitude. Models for the Fermi 130 GeV line 
must therefore explain the absence of annihilation into charged (or hadronic) modes. 

This brings us to the second interesting fact. The best fit cross-sections for the 130 GeV feature 
01E] are more than an order of magnitude smaller than the expectation for (s-wave) thermal freezeout. 
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If annihilation to SM or charged modes is absent or suppressed as the continuum limits suggest, 
reconciling this inefficient annihilation to photons with the WMAP relic density Q x h 2 = 0.1109 ± 
0.0056 [27] requires either: (1) a nonthermal relic abundance [281 129] ; (2) non-photonic annihilation 
modes which are suppressed, possibly only post-freezeout, relative to the naive radiative scaling of 
Eq. ([I]) [30, 29j; or (3) a mechanism to enhance the thermal annihilation cross-section to photonic 
final states in the early universe |31j . 

In the present work we will study two mechanisms which give enhanced DM annihilation to photons 
while also obtaining the correct thermal relic abundance, and build simple reference models for both. 
Our first example is resonant freezeout, where the presence of a resonance in the DM sector spectrum 
enhances the annihilation cross-section into photons during freezeout. Our second example introduces 
an intermediate annihilation mode, so that DM annihilation proceeds through a cascade decay of a 
non-photonic intermediate state, XX ~ > ~^ ^-1- In this example, the relation between the cross- 
section necessary for thermal freezeout and the present-day gamma line cross section is explained by 
the interplay of s-wave and p-wave contributions to the annihilation. 

In section [2] we perform an effective operator analysis of DM-photon couplings and demonstrate 
the need for new particles at mass scales comparable to the DM mass. In section [3] we perform 
a detailed examination of resonant freezeout, and consider cascade decays in section [4] Section [5] 
contains our conclusions. 

2 Effective Operators for DM Annihilation to Photons 

To demonstrate the need for multiple new states in thermal models for the Fermi 130 GeV line, we begin 
by writing effective operators to describe the interaction between DM and one or two photons. We 
assume an unbroken Z2 symmetry to explain the stability of the DM particle, and for simplicity work 
after electroweak symmetry breaking. We first consider the case when DM is a (Dirac) fermion. The 
leading operators coupling DM to photons are the electric and magnetic dipole operators, xa^xF^ 
and XPiivXF^ 1 \ appearing at dimension five [32^ I33j . However, the dominant annihilation process 
mediated by these two operators is XX ~> If through a photon in the s-channel [30]. This problematic 
annihilation to // can be suppressed in the present day if the Dirac fermion is split into two Majorana 
Xi 2j with a small mass gap m2-mi such that the depleted abundance of the heavier X2 post-freezeout 
shuts off the charged annihilation channel [301 129] ; in this case, the DM is the Majorana xi an d the 
EFT containing only xi is indeed insufficient to describe the freezeout. 

At dimension six, there are four operators which couple DM pairs to two photons, two CP- 
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conserving operators 



J^XKF^ + ^jcnxFrF*" , (2) 



and two CP-violating operators 

ij^XXF^ + ^X75Xi^ • (3) 

In addition there are operators with a single photon, such as xx^uZ^ , XlsX^nu^ 1 ^ \ an d Xl^X^^ d u h. 
We will concentrate on the operators with two photons, but we will comment on other operators when 
we introduce a concrete UV model. 

For scalar DM, the first DM-photon interactions appear at dimension 6, 

O ac , 1 = <l>%d v <l>Fi»', O sc , 2 = \<f>\ 2 F flu F^, O SCy3 = \<t>\ 2 F^F^. (4) 

For simplicity we have taken <j) complex, but this is only necessary for O ac \. As for the fermionic 
dipole operators, O sc \ will dominantly mediate annihilation to //, which can be suppressed by a 
mass splitting between real and imaginary parts of <j>. 

For fermionic DM, the presence of independent operators which contribute in different leading 
partial waves to the DM annihilation cross-section show that it is easy to accommodate an apparent 
suppression in the DM annihilation cross-section. In an EFT consisting of the CP-conserving operators 
in Eq. (j2J), the cross-section is 

4 

7YI 

crv = [4 eg + (2 4 + c\) v 2 ] + 0(v 4 ) = s+pv 2 + eV). (5) 

For a mixed partial wave freezeout process, the dark matter relic abundance is approximately given 

by m 

n h 2 „ 1-07 x 10 9 xp 

x ~ GeVM p i^f s + 3{p- s/4)/x F ' U 

in terms of the freeze-out temperature 

5 /45 g M p im x (s + 6p/xp) 



xp = In 



2vr 3 



'9 y/XF 



(7) 



In Fig. [T] we show in the left panel the region of s and p giving the correct relic density, and in the 
right panel translate that into the region of operator coefficients ci/A 3 and C5/A 3 . From Fig. [TJ we can 
see that to simultaneously accommodate the thermal relic abundance and the best-fit cross section for 
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Figure 1: Left panel: the red solid lines are the allowed region of s and p cross sections in pb to 
satisfy the dark matter relic abundance within one a. The blue dotted line denotes s = 0.042 pb-c, 
the present-day gamma line cross section in Ref. [8 J (Einasto profile). Right panel: same as the left, 
but in terms of operator coefficients ci/A 3 and C5/A 3 . The effective number of relativistic degrees of 
freedom is taken to be 75.75, g x = 4, and m x = 130 GeV. 



the Fermi-LAT 130 GeV line requires a large hierarchy between p- and s-wave scattering, p/s ~ 130, 
or C1/C5 « 23. 

Fig. [H also shows that to obtain a reasonable relic abundance through annihilation to photons only, 
the cutoffs of the effective operators are order 500 GeV if c\ and C5 are order of unity. Unfortunately 
the natural magnitude of c\ and C5 is a/n, so without additional enhancement factors, the necessary 
cutoff A is much smaller than m x , invalidating the EFT. 

In the next two sections, we study two models which extend some of the effective operators of in 
Eq. ([2]) and Eq. ([3]) by introducing a (pseudo-)scalar degree of freedom <p which couples to hypercharge 
gauge bosons through a loop of charged fermions, 

2 

C = W - 9yYB^ - m f )fi - \ r x xx<P ~ A^X75X0 - VfHifa ~ -^V , (8) 



where Y is the hypercharge of the new charged fermions fa under U(1)y , and i = 1 • • • Nf is the flavor 
index of the Nf fermions. We integrate out fa to generate the following effective operators [35] 



N f y f Y 2 a Y 2 Ai3 N f y f ae 2 f 2 



a 



- 4> D J - 4> F, V F^ = —— 4> F^F^ . (9) 

47rmj 3 47rm/3 47rm/ 

Here we have absorbed the electric charge e/, the Yukawa coupling t/f and the multiplicity factor n/ 
into the definition of There are two more operators F^ U Z^ U and Z^ U Z^ U generated. They provide 



x If 4> is a pseudo-scalar, coupling as 0i/ i 75/i, one needs to replace 2 x 4>B liu B^ v by 1 x <j)B llv B ,iv . 




Figure 2: Feynman diagram for resonant annihilation. 

a sub-leading contribution to the dark matter annihilation cross section because of the < 9^ ~ 0.05 
suppression factor. We neglect these operators in our analysis. We also point out that if the ~ 114 GeV 
gamma line becomes robust, our analysis should be extended to include operators constructed from 
W® v to fit the relative fluxes in XX ~^ 77 an d XX ~^ 

3 Resonant Freezeout 

One way to enhance the dark matter annihilation cross section at freezeout relative to the cross-section 
today is by introducing a resonance with mass slightly above twice the DM mass, 

m 2 a = 4m 2 x (l + 5) , (10) 

for < 5 <C 1 |31|, [36] . In this section we will discuss the parameter space for a resonant freezeout 
model where DM annihilation to photon pairs is entirely responsible for setting the thermal abundance. 
To be concrete, we consider a simple UV completion of the EFT in the previous section, namely Dirac 
fermionic dark matter x coupling to photons through a pseudo-scalar a with 

C D -iX\xisXa ~ JX aF ^ U ■ (11) 

We take the loop-induced aFF coupling to be given by a loop of charged fermions, as shown in Fig. [21 
so the effective coupling is 

1 = ay/JV/ej ^ _q_ 
A nm f nm / 

where yj is the Yukawa coupling of a to the heavy fermions / whose mass, number, and charge are 
given by mj, Nf, and e/. 
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3.1 s-wave 

The operators in Eq. (jlip yield the annihilation cross section 



(TV 



m x ( 


' aX x\ 
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nrhf j 


[(4m 2 + m 2 t> 2 ) — to 2 ] 2 + T 2 m 2 
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( aX x 
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1 TTlflf 


j (5 -v 2 / '4) 2 + 7 2 (l + 25) 



(13) 

where T a is the total width of a, 7 = r a /m a , and in the second line we have taken v, 6, 7 <C 1 [31 j . 
We can write the thermally averaged cross-section as 

(av) = (av) 00 f{x;5,-f), (14) 

where x = m x /T, 




1 



64vr yvrm/ ^ <5 2 + 7 2 (1 + 2<5) 

is the cross-section at x — > 00, and 



(15) 



S,) = ^r v*dv e -^/4 ^ 2 + 7 2 (l + 2^) 

n ' ' 7j yfaj (5 + ^/4)2+72(1 + 25) Ubj 

contains the information about the nontrivial velocity dependence of the cross-section. The width of 
the pseudo-scalar is bounded from below by its couplings to dark matter (in the relevant parameter 
space, the radiative width into photons is negligible). For 6 <C 1, 

We first establish that our model has a reasonable range of parameter space which can give a 
sufficiently large cross-section for XX ~^ 77- Requiring (ctv)q = (av) oo / (xo) , 

/64vr 3 m 2 \ / m f /A l v V / 1 \ 2 1 1 



I a 2 A^ i /u \ 300 GeVy VO-06/ <5 2 v ' S 2 ' 

where we have used the Einasto value for (av)o and dropped terms of order 7 2 . This determines rhf 
as a function of 5 and A^,, showing that obtaining the present-day cross-section requires 5 < 0.1. 

After using the present-day best fit cross-section to fix rhf in terms of A^, 5, the thermal freezeout 
is controlled by S and the minimum allowable width, as determined by Eq. f)lT[) . Figure [3] illustrates 
the dependence of the final relic abundance on 8 and A^,. Since the resonance sees significant overlap 
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Figure 3: The dependence of the relic abundance ^ x /i 2 on 5 for different values of X l x (left) and rhj 
(right), with the present-day cross-section fixed to the (Einasto) best-fit value for the Fermi-LAT 
gamma line excess. 



with the bulk of the velocity distribution during freezeout, the strength of the resonant enhancement 
to the cross-section is enhanced when 7 is smaller, increasing the value of the cross-section at the pole. 
Therefore if 7 is too small, e.g. A^. < 0.4 in Fig. [3J the annihilation is too efficient and the yield is too 
small to account for the present-day dark matter abundance in the absence of other decay modes for 
a. If, on the other hand, 7 is too large, e.g. A* > 0.6 in Fig. [31 annihilation is inefficient and could even 
over-close the universe. Annihilation becomes less efficient both gradually at large 5 and rapidly at 
very small 5, as the pole passes outside the main peak of the Maxwell-Boltzmann distribution during 
freezeout. 

In Fig. U] we show the contours in the S-X l x plane yielding Q x h 2 = 0.11. To illustrate the as- 
trophysical uncertainties, we quote results for both the Einasto and NFW gamma line best-fit cross 
sections. Comparing the resulting contours, we can see that the astrophysical uncertainty introduces 
an order 50% uncertainty on the couplings of the pseudo-scalar to dark matter. In the green dashed 
line of Fig. HI we show the contour for the case where a has an additional decay mode XX* with 
Br(a -> XX*) = 10Br(a -> 77). 

As can be seen from Fig. [U there is no real upper bound on the (scaled) fermion mass fhf. 
However, from Eq. (|18p. rhf is inversely proportional to S, to obtain the present-day annihilation cross 
section. If one requires 5 > 0.01 from considerations of fine-tuning, the scale for the charged fermion 
masses should be bounded from above by around 500 GeV, or in other words, within collider-accessible 
energies. 



7 




Figure 4: Contours in the ^-Al. plane (left) and rhf- X l x plane (right) yielding & x h 2 = 0.111 after fixing 
present-day cross-sections to best-fit Fermi-LAT gamma line values. The dashed line is the contour 
for the present-day Einasto gamma line cross section with the assumption that Br(a — > 77) = 10%. 



3.2 s + p-wave 

We now consider the case where more than one partial wave is important for freezeout. For concreteness 
we extend the single pseudo-scalar model of the previous subsection to include a CP-violating coupling 
X^axx to the dark matter, 

(X ~ 

C D -i\ l x xisxa ~ \ x axX ~ j-~a F^F» U . (19) 

Pure p-wave freezeout would require extremely tuned values of the resonance mass, 5 ~ 10~ 5 , to 
obtain the present-day observed cross-section. In the more interesting case of mixed s and p-wave 
resonant freezeout, the present-day observed cross-section, and hence the necessary values of 5, are 
set by the s-wave contribution, but p-wave scattering can dominate the annihilation at freezeout. 
When p-wave scattering dominates freezeout, the additional v 2 dependence of the cross-section shifts 
the velocity integral to higher values, reducing the overlap with the pole; narrower resonances are 
therefore required to achieve cross-sections sufficiently efficient to yield the observed relic abundance. 
Therefore when p-wave scattering dominates during freezeout, additional contributions to the a width 
are more constrained than in the s-wave case. Since we will be interested in cases with a hierarchy 
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between p-wave and s-wave modes, it is important to include the contribution to 7 from A£, j r = 




Figure 5: Left panel: contours in the <5-A* plane yielding Q x h 2 = 0.11, for different ratios of A£/A* . 
Results are shown for the the Einasto best fit cross-section for the gamma line excess. Right panel: 
the same as the left one but in the rhf-Xt plane. 



As for the pure s-wave case, we show the contours of f2 x /i 2 = 0.11 in the 5-Al and m/-Al planes 
in Fig. [5j We plot three different ratios of XC/XL- For A* > X^y/S, the contribution to the resonance 
width from XC, is negligible. The additional p-wave contribution to the freeze-out cross section prefers 

A 

a smaller value of A^, which explains why the pure s-wave contour (red, solid) in the left panel lies 
at larger values of A^. than the blue dashed contour, which denotes a case where p-wave scattering is 
important for annihilation but not for the resonant width. For A^, < X^^\f5, the contribution to the 
resonance width from XC, is non-negligible, reducing the enhancement from the resonance. As a result, 

A. 

larger values of \ l x are required to obtain sufficiently efficient annihilation, as can be seen from the 
green, dotted contour in the left panel. 

As can been from the right panel of Fig. [SJ the photon vertex scale rhf is again below around 
500 GeV for 5 > 0.01, implying collider-accessible charged particles. 
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Figure 6: The Feynman diagram to describe DM cascade annihilation into four photons. 

4 Cascade Annihilation 

Another way to reconcile the large photon line cross-section with the lack of continuum photons while 
obtaining the proper relic density is to have the dark matter first annihilate into (neutral) intermediate 
states, which then decay into photons. The advantage of extending the DM annihilation process with 
a cascade decay is that the relic abundance is now only controlled by the couplings of dark matter 
to the intermediate particles. The small radiative couplings of photons to the intermediate state are 
only relevant for the lifetime of that state, and are irrelevant for the annihilation cross section. The 
gamma ray signals from cascading DM annihilations have been discussed in [37|, [38l [39] and applied 
to the Fermi-LAT gamma line excess in [40|, [§]. Here, we write down a simple explicit model and 
discuss the parameter space for realizing a thermal relic abundance together with an explanation for 
the Fermi-LAT gamma line excess. 

As before, we consider Dirac dark matter x, together with a real scalar field which we denote <f>, 
and study the following set of interactions 



integrated out to generate the <pFF operator; we require only that mt is sufficiently large to forbid DM 
from annihilating into charged fermions. We are interested in the parameter space where < m x , 
so the dominant annihilation channel for DM is XX ~^ fifi ~~ ^ ^7, as shown in Fig. [6l 

The energies of the two photons from each (f> decay are, in the dark matter rest frame, 



where 9 is the angle between the photon direction in the 4> rest frame and the <f> direction of motion. 
Because 4> is a scalar field, the distribution is isotropic in 6, and the photon spectrum is evenly 





(21) 
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distributed between the kinematic endpoints: 



dN 7 
N^dRy 



which, as e 
function. 
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(22) 



0, limits to a delta function centered at m x /2. Here 0(x) is the usual Heavyside 
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Figure 7: The energy-smeared rectangular photon spectra of Eq. (|22p for different values of the mass 
splitting e. The dark matter mass is chosen to be 260 GeV for the red solid, green dotdashed, orange 
dotted lines. As a comparison, we also show a smeared delta function at 130 GeV in the blue dashed 
line. 



Since the gamma line spectrum can provide a good fit to Fermi LAT data, we also anticipate a 
good fit for this model for sufficiently small e. To estimate upper bounds for e, we consider the average 
Fermi-LAT energy resolution for gamma ray energies above 50 GeV, a(E)/E ~ 0.10 + 0.0001 E/GeV 
[41 1 . We use this energy resolution to smear the spectra in Eq. (I22p for different values of e and 
compare them with a smeared delta function spectrum centered at 130 GeV. The results are shown 
in Fig. [71 where we have shown three different values of e = 0.01, 0.02, 0.04 in the red solid, green dot- 
dashed, and orange dotted lines. The delta function spectrum is shown in the blue dashed line. As 
one can see, once C(y / e/2) is smaller than the energy resolution (order 10%) the distinctions between 
the smeared delta function and the smeared cascade spectra are minimal. 

To work out the parameter space for this model we translate the best-fit line spectra of [8] to the 
cascade model. Since the photon flux from dark matter annihilation is inversely proportional to the 



11 



square of the DM mass, and there are now four photons in the final state, the required annihilation 
cross section for this case should be approximately twice that found for a photon line. Thus we require 
{crv)o = 0.084 pb-c for the Einasto profile and 0.152 pb-c for the NFW profile. 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 8: The contours of couplings At and A^ which yield the best-fit cross-sections for the gamma 
line excess. The solid circles indicate the points which yield a thermal relic abundance of r2 x /i 2 = 0.11. 
The green dotted line is the case when Br(</> — > 77) = 0.2 of the total. 



Our explicit model considers fermionic DM and, for economy, a single real scalar eft. Thus parity 
dictates that the annihilation is p-wave suppressed unless both Ar, and the CP-violating A* are nonzero. 

A A 

In terms of the mass splitting e, the cross-section for XX ~^ fifi is 



av 



\r2\i2 



V2- 



■nrnt- 
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A^ 2 (16A 
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87A; 2 )^ 



64^' 



+ 0(v 4 



(23) 



where we have kept the leading terms in the limit e<l. 

In Fig. [HI we show the contours in the A^-A£ plane which give the required dark matter annihilation 
cross section for explaining the Fermi-LAT gamma line excess, for m x = 260 GeV and e = 0.02. The 
points which give the relic abundance ^ x /i 2 = 0.11 are indicated by heavy circles. Hierarchies of 
order ~ 10 between CP-preserving and CP-violating couplings are required, suggesting small (but not 
tiny) CP violation in a DM sector. We also show the case where the branching fraction of <f> into two 
photons is 20%. 
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An equally well-motivated alternative to the CP-violating reference model considered here would 
be to allow the Dirac DM \ to annihilate to a nearly degenerate pair of scalars with opposite parity, 
(j) and a, with subsequent decays to photons. In this case, the s-wave annihilation \X ~ * ( P a is 
proportional to (A^A^) 2 , while the p-wave annihilation is proportional to ~ (A^) 4 + {^x) 4 - If the 
branching fractions of cp, a into photons are order 1, accommodating the dark matter thermal relic 
abundance and the present-day annihilation cross section for the Fermi-LAT gamma line still requires 
a factor of ~ 10 hierarchy between A^ and A°. 

5 Discussion and Conclusions 

We have explored the possibilities for thermal models for the Fermi-LAT 130 GeV gamma line excess 
where the dominant DM annihilation channel is into photons. We consider two mechanisms, (1) 
models where a resonance in the DM spectrum enhances the annihilation rate during freezeout and to 
a lesser extent at the present day, and (2) models where DM annihilates to a new intermediate state 
which subsequently decays to photon pairs. Here the interplay of the s-wave and p-wave annihilation 
is responsible for reconciling the necessary cross-section at freezeout with the observed cross-section 
today. For the resonance model, charged fermions at the TeV scale are predicted and are accessible 
at the Large Hadron Collider, which could be interesting if the excess in the Higgs diphoton decay 
channel persists [42, 43] . 

Both of the classes of models considered in this paper require coincidences in the mass spectrum. 
For resonant freezeout, the resonance mass m a must be within a percent of twice the dark matter 
mass, a striking coincidence. For cascade decays, the intermediate state(s) <ft must be within again 
about a percent of the dark matter mass in order to have a sharp enough spectral feature to fit the 
data well. These near-degeneracies are suggestive of a dark sector with one scale A m setting the 
overall mass scale, and another Aq <C A m determining splittings. Composite dark sectors are an 
appealing avenue to flesh out the models we have considered. Near-degeneracies in the dark sector 
spectrum, as necessary for the cascade annihilation model, are readily accommodated in composite 
sectors. Resonances above threshold pose a somewhat more complicated picture. As we know from 
heavy quarkonia in the SM, a spectrum with resonances slightly above threshold could certainly be 
obtained; the complication is that such resonances would necessarily be accompanied by a resonance 
below threshold, giving the thermal relic abundance a detailed dependence on the parameters of the 
model. We leave detailed model building to future work. 
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